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Abstract. Several hybrid tomographic methods utilizing ultrasound modulation have been 
introduced lately. Success of these methods hinges on the feasibility of focusing ultrasound 
waves at an arbitrary point of interest. Such a focusing, however, is difficult to achieve in 
practice. We thus propose a way to avoid the use of focused waves through the so called syn- 
thetic focusing, i.e. by the reconstruction of the would-be response to the focused modulation 
from the measurements corresponding to realistic unfocused waves. Examples of reconstruc- 
tions from simulated data are provided. This non-technical paper describes only the general 
concept, while technical details will appear elsewhere. 

Introduction 

The last decade has seen the proliferation of the so called hybrid methods of medical imag- 
ing, where different physical types of radiation are combined into one tomographic process 
(e.g., [6, 13, 16, 26, 27, 29, 31]). Such a combination allows one to alleviate deficiencies of each 
separate type of waves, while combining their strengths. Thermoacoustic (and closely related 
photoacoustic) tomography is an example of such a hybrid method (see [2, 11, 16, 25-27, 29, 31] 
and references therein). 

In this text we discuss a very specific type of hybrid imaging, the one that combines electrical 
or optical measurements with a concurrent scanning of the object with ultrasound (US) waves. 
The propagating US wave slightly modifies the local physical properties of the medium, such 
as its electrical conductivity or the distribution of light scatterers. This, in turn, perturbs the 
background tomographic measurements (optical or electric). These perturbations are mea- 
sured and used to reconstruct the electrical or optical properties of the medium. Such hybrid 
methods hold a promise to improve the tomographic modalities that are otherwise notori- 
ous for their instability and/or low resolution, such as electrical impedance tomography [9] 
or optical imaging [29]. For these hybrid methods we will use the names Acousto-Electric 
Tomography (AET) and Ultrasound Modulated Optical Tomography (UMOT) 1 . 

It has been understood that if one could focus the ultrasound on a small spot inside the body, 
knowledge of this location would have a stabilizing effect on the reconstruction in otherwise 
highly unstable imaging modalities (see e.g. [5,7] and Section 2.2.3 below). In this text, 
however, our goal is not to discuss the ultrasound modulated imaging in detail, but rather to 
address the assumption of a well focused US beam. It is known that such focusing (to the 
extent needed for tomography) is hard to achieve (see, for example, the detailed discussion of 
this issue in [14]). Thus, it would be important to find a way of utilizing instead unfocused 
US waves. As we will show, this can be achieved through the synthetic focusing, i.e. by 
extracting the measurements corresponding to well-focused beams, from the data obtained 
with unfocused waves. Under such approach, the tomographic problem is solved in two steps 



1 Other names have also been used [6, 7, 29]. 
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(hopefully both stable): first, synthetic focusing, and then inversion from the "focused" beam 
data. We will not attempt to address the second step in this short text (see [5, 7, 17] where the 
various implementations of such a procedure are discussed). Instead, we show that the first 
step (synthetic focusing), under appropriate choice of waves, is stable and mathematically 
equivalent to the reconstruction in thermoacoustic tomography (although no thermoacoustic 
measurements are conducted). In some cases, this step can even be reduced to the inversion 
of the Fourier transform. An example of the AET reconstruction is also provided. 

This paper is meant as a preliminary non-technical announcement, which follows a part 
of the lecture given at the Jan Boman's conference "Integral Geometry and Tomography" in 
August 2008. The technical analytic and numerical details will be provided elsewhere. 

The structure of the paper is as follows: The idea of US modulation in tomography is 
summarized in Section 1. Section 2 contains the description of four different types of unfocused 
US waves that one might try to use. In each case it is shown that the synthetic focusing (within 
the assumed mathematical model approximation) reduces to one of the well known analytic 
reconstruction procedures. Section 2.2.3 is devoted to a brief description of the AET procedure 
and a numerical example. Remarks are contained in Section 3, followed by Acknowledgments 
and References. 

1. Ultrasound modulation in tomography 

As indicated in the Introduction, the ultrasound (US) modulation in various types of to- 
mography is performed by sending a US wave through the object, concurrently with the 
original tomographic measurement (say, optical or electric). Varying the shape of the wave, 
one observes the resulting changes in the tomographic data and tries to extract from them 
the quantity of interest. This technique promises to improve significantly the stability of such 
modalities as electrical impedance tomography (EIT) and optical tomography (OT). 

To be more specific, let us consider the EIT case 2 , where one seeks to recover the conductiv- 
ity a(x) inside of an object occupying the domain Q from boundary impedance measurements. 
Assume for instance that a boundary current g(y) is applied and one measures the correspond- 
ing boundary potential response h(z). Here y and z denote variable points on the boundary 
dQ. In EIT, the current g is varied and the potential h is measured, so that the complete 
Neumann-to-Dirichlet operator A CT on dQ is obtained (alternatively, by varying h, one recovers 
the Dirichlet-to- Neumann operator). In the case of ultrasound modulation (at least in the 
example discussed in Section 2.2.3), a single boundary current can be used, rather than the 
whole operator A CT ; we will thus assume that g is fixed. Suppose now that, given any point 
x G f2, we could create a US wave that would approximate well the delta function at the 
location x. This would create a perturbation h x (y) of the original boundary voltage h(y). By 
scanning the focusing point x through the domain Q, one obtains the set of functions h x , x G Q, 
from which one can try to recover the internal conductivity a. This can actually be done, 
and the reconstruction does not inherit the original high instability of EIT (see the details in 
[7, 17], where two different approaches are presented). A somewhat similar situation arises in 
OT [5,23,24]. Stability of such reconstructions can be seen in Figures 1 and 2 (taken from 
[5, 17]). These figures show mathematical phantoms (left) and their reconstructions obtained 
using the simulated perfectly focused US modulations in EIT and OT correspondingly (images 
on the right). These examples clearly demonstrate stability of these hybrid modalities; the 
issues of stability will be discussed in more detail, along with the reconstruction methods, in 
[5,17]. 



2 The first AET reconstructions were probably suggested and implemented in [32]. 
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Figure 1. Phantom (left) and its EIT reconstruction (right) using focused US modulation. 




Figure 2. Phantom (left) and its OT reconstruction (right) using focused US modulation. 

However, as we have mentioned in the Introduction, the problem with this technique is the 
difficulty of physically generating well localized US waves (see e.g. [14]). We thus suggest the 
use of a synthetic focusing approach described in the next section. 

2. Synthetic focusing 

We start by mentioning that the acousto-electric and opto-acoustic effects that are respon- 
sible for the perturbations in the electric and optical properties of the medium due to US 
irradiation, are very small (e.g., [19,20,29,32]). Although this makes the resulting boundary 
effects harder to measure, it also enables one to use the linearization of the problem in the 
following sense. Let a(x) be the parameter to be reconstructed, say, electrical conductivity 
in the EIT case. The perturbation 5a of a due to the applied US wave will be very small. 
Hence, one can safely assume that the operator L that maps the perturbation 5a into the 
perturbation 5h of the measured data h is linear. In fact, the operator L depends on the 
underlying a and thus will be denoted as L a . 

Due to smallness of the amplitude of the applied US wave, as well as the smallness of 
the acousto-electric and opto-acoustic effects, one can try also to linearize the effect of local 
pressure on the perturbation 5a. It turns out that such a linearization is possible in AET; it 
happens to be incorrect in the optical case, where the dependence is quadratic [29]. 
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Let now 

(1) L a : 5a(x) i-> 5h(y), x efl,yeTcdfl 

be the linear operator that maps the distributed perturbation 5a{x) of the material parameter 
(say, conductivity) inside fl to the perturbation 5h(y) of the tomographic data measured on 
a part T of the boundary dfl (which, in principle, could consist of a single point, although 
this would not be advisable in tomography). If one could create US waves well localized at 
arbitrary locations x G fl, this would allow one to measure the response of L to delta-type 
inputs, and thus would lead to a direct measurement of the kernel l(x, y) of the operator L: 

(2) (Lf)(y) = J l(x, y)f(x)dx, x G fl, y G V. 

n 

So, the use of perfectly focused US waves is equivalent to measuring the kernel l(x,y). The 
general idea of synthetic focusing is to use a more realistic non- localized complete set w a (x) 
of US waves and to recover the kernel l(x, y) from the measured responses 

l(x, y)w a (x)dx. 



n 



In other words, one reconstructs numerically the would-be response of the tomographic mea- 
surement to the perfectly focused US irradiation, which explains our use of the name "synthetic 
focusing" . 

We now address some specific versions of synthetic focusing in sub-sections below. (The 
reader should be warned that we discuss here only the reconstruction of the kernel l(x, y), but 
not the procedure of the further recovery of a(x) from the kernel.) 

2.1. Spherical waves. Suppose that one places a point US transducer at arbitrary location 
z G dfl and creates a short impulse that approximates a spreading spherical wave p(t,z,x) = 
S(\x — z\ — vt), where v is the constant US speed in the tissue 3 . Then the response of the 
measurement operator L a produces the spherical integrals of the kernel l(x,y) 

(3) L a (p(t, z, •)) = J l(x,y)dx,z G dfl. 

\x— z\=vt 

We thus arrive to the problem of recovering a function l(x, y) of x (y is now just a fixed param- 
eter) inside domain fl from known integrals from / over all spheres centered at the boundary 
dfl. This problem has been studied extensively, in particular due to its relevance for the 
thermoacoustic tomography (TAT); numerous results concerning the solvability, stability, and 
inversion methods have been obtained (see [1-4,10-12,16,18,25,28,29,31]). In particular, 
when dfl is a sphere, a variety of backprojection type formulas have been derived. We provide 
below the 3D formula derived in [30] (it was generalized to all dimensions in [18]): 



(4) f(x) = ±-dW [n( Z ^ ld9 ^ t) 



8vr 2 J \ t dt 



dA(z). 

t=\y-x\ 



Here fl is the unit sphere, g(y, t) is the average of / over the sphere centered at y G dfl and of 
radius t, n(y) is the outward normal vector to dfl at y, and dA is the surface area measure. 
Several different inversion formulas were derived earlier in [10]. 



3 This certainly requires a broadband transducers, as is the case in all other suggested synthetic focusing 
methods. 
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In the case of non-spherical surfaces, eigenfunction expansions and time reversal methods 
can be used to reconstruct I (see [1,2,4, 10-12] and references therein). 

2.2. Spherical monochromatic waves and plane waves in AET. Since in AET the 
electric response of the medium is proportional to the perturbation in the pressure, one can 
use for scanning not only acoustic pulses, but also monochromatic waves. Using such waves 
has the advantage that the resulting boundary measurement will oscillate with ultrasound 
frequency. Hence, its time Fourier transform allows one to selectively pick up this frequency, 
thus reducing the effects of the wide spectrum noise. 

2.2.1. Spherical monochromatic waves. Suppose that the monochromatic wave p\(t, z, x) with 
a time frequency A is produced by a point transducer located at the point z G dfl: 

e i\\x-z\ 

Px(t,z,x) = e~ iv 



Ait\x — z\ 
Then 

/g«A|x — 2 
l(x, y)e- iXt — rdx, x G Q, y G T, z G dfl, 
4n\x — z\ 

Q 

and the Fourier transform with respect to time Lp x (z,y) of Lp\ will be equal to the integral 
of l[x, y) multiplied by the free-space Green's function $\(x, z) of the Helmholtz equation: 



Lpx(z,y) = J l(x,y)$ x (x,z)dx, 



n 

p iX\x- 

<f>x{x,z) = 



A-k\x — z\ 

As before, we assume that the point y at which the electrical measurements are conducted, 
is fixed. Then, if the measurements are repeated for a wide range of frequencies A and for 
all z G dfl, and are followed by the inverse Fourier transform (in A) applied to Lp\{z,y), 
one recovers the integral of / over the spheres centered at z as in (3). Now the problem of 
recovering l(x, y) from values of Lp x (z, y) is equivalent to inverting its spherical mean Radon 
transform of /; such an inversion was discussed in the previous section. 

Alternatively, since the first step of the reconstruction formulas presented in [18] consists in 
computing the values of Lp x (z, y), one can reconstruct l(x, y) without explicitly reconstructing 
the values of its spherical mean Radon transform first. In particular, by using the 3-D formula 
from [18], one obtains the following reconstruction formula for l(x,y) from Lp\(z,y): 

l(x,y) = -^div x J n(z)h y (z,\x - z\)dA(z), 

\z\=R 

with 

h y (z, t) = —j J cos(At) Im ^Lp\(z, y)j — sin(At) Re (^Lp\(z, y)j XdX. 

Here, as before, n(y) is the outward normal vector to dfl at y, and dA is the surface area 
measure. 
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Figure 3. AET reconstruction using spherical monochromatic acoustic waves. 

2.2.2. Plane waves. Spherical waves arise when the size of the transducer is small compared 
to the wavelength. If the distance from the transducer to the object of interest is much larger 
than the size of the object, the wave is approximately planar. Plane waves can also be created 
by using a large planar transducer. In either case the transducer should be broadband to 
permit generation of plane waves in a wide range of frequencies. 

Measurement done using planar acoustic waves Pk(t, z, x) = e~*' fc '* exp(i/c • x) correspond to 
measuring the Fourier transform of the kernel l(x, y): 

Lp k (z, y) = J exp(ik ■ x)l(x, y)dx. 
n 

Synthetic focusing now reduces to the inversion of the Fourier transform. 

2.2.3. Example of AET. Let us illustrate by a numerical example the application of synthetic 
focusing to image reconstruction in AET. We simulated numerically a square domain, with 
the electrical currents equal 1 on the left and right sides of the square and on the top and 
the bottom. The conductivity a(x) we used in our experiments was close to 1; the gray scale 
density plot of the logarithm log cr(x) is shown in the left part of Figure 1. In this figure 
the light circles correspond to the value logcr(x) = 0.05, the dark ones represent the value 
of log cr(x) = —0.05, and the gray background depicts \oga(x) = 0. The simulated electric 
potential was "measured" on the whole boundary of the square. 

As mentioned previously, Figure 1 (right) presents the result of AET reconstruction from 
simulated measurements corresponding to perfectly focused US modulation. We compare this 
to the results of AET reconstruction of the same phantom, obtained by using synthetic mea- 
surements, as shown in Figure 3. In this example we modeled the perturbations of electric 
potential on the boundary caused by the spherical pulse waves with the centers on a circle sur- 
rounding our square domain. The total of 300 transducers and 800 different radii of outgoing 
spherical pulses per transducer were simulated. The reconstruction shown in the right half of 
Figure 3 was obtained by using synthetic focusing to obtain l(x, y) and then by reconstructing 
conductivity o~(x) (again, the methods for reconstructing a from I will be discussed elsewhere.) 
The reconstructed image is clearly as good as the one obtained using the perfectly focused US 
modulation (Figure 1, right part). 

2.3. Elongated focusing areas and X-ray transform. One of the UMOT reconstruction 
methods, suggested in [22] (see also [29]), is in fact a synthetic focusing, rather than a true 
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reconstruction. The authors of [22] employed the observation that the focusing was non- 
perfect: the focal zone of the ultrasonic waves was 2 mm across and 20 mm in length, which 
made it a better approximation to a segment rather than a point. One thus can conclude that 
by shifting and rotating the focusing area, one can extract the X-ray transform of the kernel 
l(x,y) (as before, with respect to the variable x). Then standard X-ray inversion formulas 
recover the kernel /. 

3. Remarks 

(1) As it has already been mentioned, exact analytic and numerical details of the recon- 
struction procedures will be provided elsewhere. 

(2) Practical implementation of the suggested synthetic focusing procedures will have to 
include more precise information about possibilities and imperfections of the currently 
available transducers. 

(3) The UMOT situation differs from AET and will probably not allow usage of monochro- 
matic US modulation. However, the modulation by spherical pulse waves should still 
be possible. These issues will be addressed elsewhere. (As it has been already men- 
tioned, the synthetic focusing using X-ray transforms was implemented for UMOT in 
[22,29].) 
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